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ABSTRACT 

Surface longwave radiation fluxes dominate the energy budget of nighttime polar regions, yet little is known 
about the relative accuracy of existing satellite-based techniques to estimate this parameter. We compare eight 
methods to estimate the downwelling longwave radiation flux and to validate their performance with measure- 
ments from two field programs in the Arctic: the Coordinated Eastern Arctic Experiment (CEAREX) conducted 
in the Barents Sea curing the autumn and winter of 1988. and the Lead Experiment performed in the Beaufort 
Sea in the spring of 1992. Five of the eight methods were developed for satellite-derived quantities, and three 
are simple parameterization* based on surface observations. All of the algorithms require information about 
cloud fraction, which is provided from the NASA-NOAA Television and Infrared Observation Satellite (TIROS) 
Operational Vertical Sounder (TOVS) polar pathfinder dataset (Path-P); some techniques ingest temperature and 
moisture profiles (al .o from Path-P); one-half of the methods assume that clouds are opaque and have a constant 
geometric thickness of 50 hPa. and three include no thickness information whatsoever With a somewhat limited 
validation dataset, tie following primary conclusions result: 1) all methods exhibit approximately the same 
correlations with mt asurements and rms differences, hut the biases range from -34 Wm • (16# of the mean) 
to nearly 0: 2) the * rror analysis described here indicates that the assumption of a 50-hPa cloud thickness is 
too thin by a factor of 2 on average in polar nighttime conditions: 3) cloud-overlap techniques, which effectively 
increase mean cloud thickness, significantly improve the results: 4) simple Arctic -specific parametcrizations 
performed poorly, frobably because they were developed with surface-observed cloud fractions whereas the 
tests discussed here used satellite- den ved effective cloud fractions: and 5) the single algorithm that includes an 
estimate of cloud th ckness exhibits the smallest differences from observations. 


1. Introduction 

In this investigation we extend the work of Key et 
al. (1996) and evaluate the ability of several existing 
methods to estimate the surface downwelling flux of 
longwave radiation (DLF) over snow- and ice -covered 
surfaces, particularly at night. The algorithms examined 
by Key et al. (1996) are al simple parametcrizations 
empirically derived from si rface measurements in the 
Arctic; here we also evalu ite several techniques that 
ingest retrievals from satell te data. Our goal is to de- 
termine which method(s) should be used in computing 
DLF from satellite retrieva s for a variety of applica- 
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lions, including a long-term dataset for studies of var- 
iability and change and use as forcing fields for mod- 
eling studies. We also extend the work of Key et al. 
(1997) in evaluating the sensitivity of DLF to pertur- 
bations in various atmospheric parameters with the in- 
tent of identifying variables that require improved re- 
trieval accuracy. 

Polar regions are notoriously problematic for global 
surface radiation budget (SRB) algorithms based on sat- 
ell ile data, because several factors complicate the re- 
trieval process and few validation data are available. 
Because the high latitudes are recognized as climatically 
sensitive areas, there is a strong demand from the sci- 
entific community for reliable, long-term surface radi- 
ation datasets for the polar regions. We employ a com- 
bination of satellite retrievals from the National Aero- 
nautics and Space Administration (NASA)-National 
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Oceanographic and Atmospieric Administration 
(NOAA) Television and InfrareJ Observation Satellite 
Operational Vertical Sounder (TOVS) polar pathfinder 
(hereinafter, Path-P) dataset (F rancis and Schwciger 
2000) with surface observations from radiometers, ceil- 
ometers, and surface observers to evaluate the perfor- 
mance of several existing algorithms to estimate surface 
longwave fluxes in the Arctic. 

Longwave radiation dominates the Arctic surface en- 
ergy budget for almost one-half of the year when in- 
solation is absent or weak. During the polar w inter, tur- 
bulent fluxes are small in sea-ice -covered regions, ex- 
cept over cracks in the ice when vertical air-ocean tem- 
perature and moisture gradients are large, fn contrast to 
lower latitudes, at which low-lex el temperature and wa- 
ter vapor content largely govern the DLF, clouds play 
the most important role in polar regions. Sensitivity 
studies by Key et al. (1997), Fiouin et al. (1988), and 
Chiacchio (2001) show that DIP is most sensitive to 
cloud fraction, cloud thickness (or liquid waiter path), 
and cloud-base height. Of these parameters, passive sat- 
ellite sensors can be used to estimate only cloud fraction 
and cloud-top height during polar night conditions, and 
even these have much larger uncertainties than do es- 
timates from other parts of the globe. Algorithms to 
detect clouds and to diagnose their properties often fail 
over snow- and ice-covered areas because of frequent 
surface-based temperature inversions that confound sat- 
ellite eloud-detcction algorithms and introduce uncer- 
tainty into satellite-retrieved temperature profiles. Short- 
wave channels, especially new' sensors that measure ra- 
diances in the 1. 6- /urn wavelength region, add consid- 
erable information, but historica visible data are limited 
by the lack of contrast between clouds and snow. Efforts 
are underway to infer polar cloud characteristics beyond 
fraction and cloud-top height, Fut they are still exper- 
imental. 

Several algorithms and parameterization* exist for es- 
timating DLF from satellite-derived information, hut in- 
tercomparisons and validation for polar-night conditions 
have not been performed. The objective of this inves- 
tigation is to compare DLF values computed with eight 
different methods quantitatively to validate results with 
measurements from surface-bas;cf instruments and hu- 
man observers, to identify prob;.ble causes for errors in 
each method, and to make recommendations as to which 
algorithm(s) provides the best estimates of DLF in the 
Arctic night. 

2. Data sources and tools 

a. Satellite 'derived product s 

Several of the methods undtr investigation require 
temperature protiles, humidity profiles, surface temper- 
ature, cloud fraction, and/or cloud height. For this study, 
atmospheric state information is obtained from the 
NASA-NOAA TOVS Path-P dataset (Francis and 


Schw'eiger 2000). The TOVS instrument, which has 
flown on NOAA polar-orbiting sensors since 1978. com- 
prises the High-Resolution Infrared Radiation Sounder 
(HIRS). the Microwave Sounding Unit (MSU), and the 
Stratospheric Sounding Unit (SSU). Data from SSU are 
not used to create the Path-P dataset. HIRS measures 
radiances in 20 channels from the visible to infrared 
wavelengths with a resolution of 17 km at nadir, and 
MSU has four channels in the oxygen absorption band 
near 50 GHz. The Palh-P dataset was produced using 
the improved initialization inversion (31) processing al- 
gorithm for TOVS radiances, developed by the Atmo- 
spheric Radiation Analysis Group at the Laboratoire de 
Meteorologie Dynamique in Palaiseau, France (Che'din 
ct al. 1985). The 31 algorithm was modified to improve 
retrieval accuracy over sea ice and snow [Francis 
(1994); Scott et al. (1999)]. Path-P products are pro- 
vided daily on a ( 100 km) : grid and include temperature 
and moisture profiles, surface skin temperature, cloud 
fraction and height, and a variety of other parameters. 
The 31 algorithm has at its core a comprehensive library 
of global atmospheric profiles (-1800) that provides 
the first guess to this physical-statistical technique and 
consequently is able to capture the strong surface-based 
and elevated temperature inversions that are nearly ubiq- 
uitous in all seasons but summer in the Arctic region. 
Validation of surface and 900-hPa temperatures with 
radiosonde data reveals small mean errors (1.4 and 2.5 
K). Retrieved inversion strength, however, is often less 
than radiosonde values owing to the coarse vertical res- 
olution of the temperature profile, and the cap may be 
misplaced in the vertical by a few tens of hectopascals. 
The cloud fraction variable (labeled FCLD in the .Path- 
P dataset) is an effective cloud fraction A.z?. which is 
the product of the fraction A, of the sky covered by 
cloud and the cloud emissivity e. Cloud emissivity rang- 
es between 0 an I ; therefore A, s is always less than or 
equal to A, , which is the quantity reported by human 
observers. This distinction is significant in polar regions 
because optically thin clouds — even in the infrared — 
are common, especially in winter. Hereinafter we ab- 
breviate effective cloud fraction as A et . and “cloud frac- 
tion" denotes the fraction of the sky covered by cloud 
(A, ). See Schw'eiger et al. (2001) for additional infor- 
mation and validation results for the Puth-P data. 

b. Validation datasets 

A significant problem in studying cloud or surface 
characteristics in this region is (he paucity of measure- 
ments, especially in winter when DLF is the dominant 
component of the surface energy budget. Observations 
from two field experiments are used in this study. The 
first is the drift phase of the Coordinated Eastern Arctic 
Experiment {CEAREX; CEAREX Drift Group 1990), 
which was conducted in the eastern Arctic Basin from 
September of 1988 through January of 1989 (Fig. I). 
The experiment included two research vessels and an 
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array of surface camps at which a variety of meteoro- 
logical and oceanographic measurements were made. In 
this study we use data obtained aboard the R/V Po- 
larbjorn , which include downward infrared fluxes (be- 
tween 4 and 50 /un) from an Eppley Laboratory, Inc., 
pyrgeometer, which has a nominal instrument error of 
5 W m 2 (CEAREX Drift Group 1990). Owing to the 
lack of solar radiation during CEAREX, as well as the 
low sun angles and large cloud fractions during the 1992 
Beaufort Sea Lead Experiment (LeadEx). we assume 
errors in radiometer measurements resulting from solar 
contamination are small. The radiometer domes required 
cleaning hourly to remove frost and precipitation; only 
measurements following a cleaning were used to com- 
pute daily-average flux values, which we compare with 
daily-average satellite-derived DLFs. The differing 
space scales of Path-P data and surface point measure- 
ments is a possible source of error. Schweiger et al. 
(2001) analyzed correlations between time-averaged 
point measurements of cloud fraction and spatially av- 
eraged satellite values and found that correlations were 
low for timescales shorter than 2 days and peak at 8 
days. Because clouds are the dominant factor in deter- 
mining DLF in the Arctic winter, we assume these cor- 
relations also apply to DLF. They speculate that the lack 
of strong correlation at short timescales may be caused 
by the differing perspectives of satellite versus surface 
observations (view from above or below). Another prob- 
able cause is that smaller variations occur at short time- 
scales, which may cause this signal to be lost in noise, 
whereas large variations may occur at long timescales 
and so are more detectable above the noise. 

Data from LeadEx (LeadEx Group 1993) are also 
used to validate DLF computed from each of the eight 
algorithms. This held program was conducted in the 
Beaufort Sea at a camp on the pack ice that drilled 


westward from 24 March 1992 until 25 April 1992 (Fig. 
I ). The main objectives of the experiment were to study 
the cracklike openings (leads) in sea ice formed by the 
ice deformation and to understand the effects of leads 
on the polar ocean and atmosphere. In addition to ra- 
diation measurements, we use cloud-base height retriev- 
als from a lidar ceilometer to compare with Path-P- 
derived values. The vertical resolution of the ceilometer 
retrievals was 30 m, and the instrument could observe 
cloud bases up to 8 km. These ceilometer estimates are 
not considered to be reliable for absolute validation, 
however, owing to reported problems in detecting thin, 
low-level, ice clouds (O. Persson 1999, personal com- 
munication). Thirty-second ceilometer values are av- 
eraged over 24 h to be consistent with Path-P products. 

c. Radiative transfer model 

Sensitivity tests and calculations of surface fluxes are 
performed using a forward radiative transfer model 
called Streamer, which was assembled by J. Key and A. 
Schweiger (Key and Schweiger 1998). Streamer is a 
publicly available, highly flexible package that can be 
used to calculate shortwave and longwave radiances and 
fluxes for a wide variety of atmospheric and surface 
conditions. Absorption and scattering by gases is pa- 
rameterized for 24 shortwave and 105 longwave bands. 
Built-in data include water and ice cloud optical prop- 
erties, aerosol profiles, and seven standard-atmosphere 
profiles, or users can provide their own. Each scene can 
include up to 10 cloud types, up to 10 overlapping cloud 
sets of up to !0 clouds each, and up to three surface 
types. Also, spectral albedos for various surface con- 
ditions are included. The number of streams used in the 
calculation can be varied; two are used in this study. 
Modeled fluxes for standard atmospheres were com- 
pared with calculations by approximately 37 other mod- 
els presented in the report of the Intercomparison of 
Radiation Codes in Climate Models project (Ellingson 
et al. 1991 ). Streamer-computed fluxes were within 57c 
( 1 standard deviation) of the mean of all models (Francis 
1997). 

3. Sensitivity of DLF to atmospheric parameters 

Sensitivity studies are performed to determine the 
likely errors in DLF resulting from uncertainty in cloud 
parameters and from published uncertainties in Path-P 
products. DLF is calculated using Streamer for typical 
winter and summer Arctic conditions and with expected 
rms errors (in parentheses) for each of (he following 
slate variables under clear conditions: surface skin tem- 
perature (±5 K). temperature profile (±3 K at all lev- 
els), and moisture profile (±30% at all levels). In ad- 
dition, we estimate the sensitivity of DLF to varying 
hulk cloud properties: fraction, thickness, base height, 
liquid water content (LWC), and effective droplet radius 
r . Each calculation includes o/one amounts for a stan- 
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Pvrturbotton to Sfc. r«mp«rotur* (K) 



Porturbotlon to Tonip. Prof ho (K) 

Fir*. 2. Sensitivity of DLF to perturbations in the (a} water vapor 
profile (30^1. (h) surface temperature 3 Kt. and (c J temperature 
profile (3 K> lor standard Arctic sumi ier (solid line) and winter 
(dashed) conditions. Dotted lines indicate 10 W m ’ error bounds. 


Tabi i I. Resulls from sensitivity tests. 


Variable name 
(ex peeled rms error in 
Path-P retrieval) 

DLF error 
(W m -> 

DLF error 
( W m ') 

Surface skin temperature |3 K) 

4 t winter 

8 (summer) 

Temperature profile (3 K) 

5 t winter) 

10 (summer) 

Water vapor profile (300* ) 

4.5 (winter) 

10 (summer) 

Cloud fraction (3fK* ) 

14 ilow cloud) 

X (high cloud) 


(Jard subarctic winter profile, and the carbon dioxide 
concentration is fixed at 350 ppmv. Clouds in these tests 
are composed of spherical water droplets with a nominal 
cloud thickness of 50 hPa, an LWC of 0.20 gm and 
a typical r, of 8 ^im (Curry and Ebert 1992). We feel 
justified in considering only water clouds, because liq- 
uid water was detected in over one-half of the clouds 
during the Surface Heat Budget of the Arctic (SHEBA) 
experiment in every month except December (Intrieri 
et al. 2001) and because phase alone has a negligible 
effect on DLF (Francis 1999). 

Results of sensitivity tests for temperatures and water 
vapor are shown in Fig. 2 and are summarized in Table 
I. These results are consistent with those of Key et al. 
( 1997) and show that errors in DLF arising from doc- 
umented uncertainties in satellite-derived temperature 
and moisture profiles will be well within the 10 W m : 
threshold that has been suggested by the World Climate 
Research Programme (WCRP) as the target accuracy 
for surface flux estimates (Raschke et al. 1990). We also 
use Streamer to test the sensitivity of DLF to uncer- 
tainties in TOVS-retrieved surface-based inversions. W'e 
calculate DLF for a typical winter surface-based inver- 
sion (13 K difference between the surface and top of 
the inversion at 900 hPa) and for a temperature profile 
with the inversion smeared out and a positive lapse rate 
throughout. A typical 50-hPa-thick water cloud is placed 
with its top at 900 hPa (top of the inversion). The dif- 
ference in DLF between these two model runs is less 
than 3 W m - (inversion run is smaller). Because this 
scenario is likely a worst case, we conclude that any 
errors in TOVS-retrieved inversion strength, height, or 
even existence would not result in the magnitude of DLF 
deficiencies we observe in many of the algorithms. 

Figure 3 shows the computed sensitivity of DLF to 
cloud fraction and geometric thickness. A 3 (Wr error in 
low-cloud fraction would result in DLF errors in excess 
of the WCRP threshold, but DLF is less sensitive to 
uncertainties in high-cloud fractions. DLF is highly sen- 
sitive to cloud thickness — more so for low clouds: an 
error of 10 W m ; could arise from assuming a cloud 
is only 20 hPa thicker than its actual value. This is an 
important point in our later discussion of the assumption 
in some algorithms that clouds are a constant 50 hPa 
thick. 

Figure 4 shows the sensitivity of DLF to cloud LWC 
for varying droplet effective radii and cloud-base 
heights, with and without a surface -based temperature 
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Cloud thickness [mb] 

l*Ki. 3. Sensitivity oT DLF to perturbations in (a) cloud fraction 
and (b) cloud thickness, with cloud-base heights at 800 and 400 hPa 
in typical Apr conditions in the central Arctic. Clouds are composed 
of spherical water droplets w ith a nominal cloud thickness of 50 hPa. 
an LWC of 0.20 g m *. and a typical equivalent radius r of 8 /im. 

inversion. We use temperature and water vapor profiles 
typical for April, and LWC is varied from 0.0 to 0.2 g 
m \ In each of these experiments, it is apparent that 
the DLF is extremely sensitive to the LWC in thin, high 
clouds that contain less than 0.02 g m 1 of water and 
in low clouds with less than 0.05 g m \ For thicker 
clouds, DLF is no longer sensitive to LWC; that is, the 
cloud is optically thick in the infrared. Cloud droplet 
effective radius in this size range, however, has a neg- 
ligible effect oil DLF. This result is consistent with re- 
sults by Francis (1999) that show little sensitivity of 
infrared cloud radiative forcing by water clouds (r r = 
10 /xm) versus ice clouds (r^ — 50 gm). Figure 4b shows 
that DLF is more sensitive to the LWC of low clouds 
than of high clouds and that, for optically thick clouds, 
an error of 10 W m " 2 would result from an error in 
cloud-base height of approximately 150 hPa. Figure 4c 
illustrates the effect of a cloud base lying above and 



0,00 0.05 0.10 0.15 

LWC [gm"*] 

Fi(i. 4. Sensitivity of DLF to perturbations in LWC with (a) r. = 
4 and 10 /am, (b) cloud-base heights at 400 and 800 hPa with no 
temperature inversion, and (c) eloud-base heights at 860 and 950 hPa 
with a temperature inversion. Ordinate is cloudy minus clear DLL 

below the cap of a surface-based temperature inversion. 
The results of these sensitivity tests for bulk cloud pa- 
rameters are also generally consistent with those of Key 
et al. (1997), although we test some different variables 
and, in some cases, over a more widely varying range 
of values. We are interested only in the sensitivity of 
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I \hi i 2. Summary ol inputs ami assumptions tor Dl.F algorithms. 


Algorithm reference 

Required input data* 


Assumptions* 


Francis ( 1997) 

</(;>. A f\ TBs. T. 

One cloud layer 



Frouin et al. (19X8). FI with overlap 

Ti:). </(;». A /’, / , 

Overlapping cloud layers. / - 50 hPa. r. 

- 1. Ttz) < 0 

Frouin et al. (1988). F2 

/ 1 : ). qi : ). A 

One cloud layer, 7 

- 50 hPa. ft — 1 . 7CI 

l < 0 

Marshunova (!%6) 

A. T <' 


— 


Mavkul and Church (1973) 

A, T. 


.... 


Gupta (19891 

7*1 ;). i/C). A t\ 7\ 

One cloud layer. Z 

= 50 hPa. t: = l./<;l 

<. 0 

Zillman (1972) 

A f\ 


— 


Frouin et al. ( 1988), FI 

/ U-). qi:), A r 7\ 

One cloud layer. 7. t 

50 hPa. k = 1. 7 ( r 

1 < 0 


* Here 7‘(-) and qi-) arc temperature and moisture pr- dries, A, is cloud traction. A, is effective cloud traction. / > , is cloud-top pressure. 
Z, is cloud thickness, t: is cloud emiss vily. /'. is surface temperature, v is near-surlace vapor pressure, and TBs are HIRS brightness 
temperatures. 


DLF to individual variables, given that part of our goal 
is to identify which satellite-retrieved quantities lack 
sufficient accuracy for DLF ealeu at ions. For an analysis 
of the overall uncertainty in DLF see Key et al. (1997). 

From the results of these tests, we conclude that DLF 
is most sensitive to errors in clou j fraction and to LWC 
in thin clouds and that DLF' is relatively insensitive to 
droplet size. Known uncertainties in satellite-retrieved 
temperature and water vapor profiles result in DLF er- 
rors within 10 W m % For further details on sensitivity 
tests, see Chiacchio (2001). 

4. Methodology 

We evaluate the ability of seveial methods to estimate 
DLF by comparing daily-mean-calculated values from 
a (100 km)- grid box to surface -measured DLF from 
two held experiments in the Arctic. Some of the algo- 
rithms, such as those of Gupta (1989) and Frouin ct al. 
(1988), are being used globally, and we compare their 
performance with Arctic-specifi : algorithms, such as 
that of Francis (1997). We also evaluate throe simple 
empirical parameterizations developed for Arctic con- 
ditions: Marshunova (1966), /ill nan ( 1972), and May- 
kut and Church (1973). These parameterizations are 
among the longwave flux methods examined by Key et 
al. (1996). yet we do not consider the Schmelz et al. 
(1986) method, because it was developed for daylight 
conditions. In this section we describe each method and 
the required input data. A summary of the assumptions 
and required information for each method is given in 
Table 2. 

a. Gupta (1989) 

This algorithm was used to generate a global 8-yr 
SRB dataset (Gupta ct al. 1999* and also is included 
among other algorithms in both the Clouds and the 
Earth's Radiant Energy System (CERES) and the 
WCRP-GIobal Energy and Wa er Cycle Experiment 
(GEWEX) SRB projects. This method requires inputs 
for water vapor and temperature profiles, cloud fraction, 
and cloud-top pressure, which Gupta (1989) obtained 
from NOAA operational TOYS retrievals, and the 


GEWEX SRB and CERES projects obtained from re- 
analysis datasets (either the European Centre for Me- 
dium-Range Weather Forecasts or NASA's Data Assim- 
ilation Office). In this study, we instead use the Path-P 
products, because they are believed to be more accurate 
in Arctic conditions (Francis 1994). The primary as- 
sumptions of the Gupta ( 1989) technique are that clouds 
exist in a single. 5()-h Pa- thick layer and they are opaque. 
The method parameterizes DLF as 

DLF - DLF ,(1 - A ) + DLF c)llv A,, (1) 

where DLF fc . lf is the downward longwave flux for clear 
sky, DLF ck)v is the downward flux for cloudy sky, and 
A, is the cloud fraction, for which we use A^ from Path- 
P. This equation is then simplified as 

DLF = C, + C\A IV , (2) 

where C, = DLF,, and C : = (DLF Mn - DLF U , ). that 
is, the thermal emission from the cloud base. The C, is 
parameterized in terms of the effective emitting tem- 
perature F of the atmosphere (estimated front the sur- 
face and lower tropospheric temperatures) and the sur- 
face-to-700-hPa water vapor burden W, (mm): 

c, - f{W x ) 7';., (3) 

where x is empirically determined to be 3.7. Further. 
/CIV,) is expressed as 

/( VV, ) - C,/F; - A n + A,V + A 2 V 2 + A,V\ (4) 

where V = ln( VV', ) and the As are regression coefficients. 
Equation (4) is then fit to fluxes and meteorological 
profiles from five sites in the United States to obtain 
the regression coefficients A 0 = 1.79 X 10 A, = 
2.093 X 10 LA, = -2.748 X 10 Land A, = 1.184 
X 10 L 

The parameterization of CL contains the terms F.,, 
(cloud-base temperature) and IV (water vapor below the 
cloud). To determine F h , the cloud-base height must be 
calculated from the satellite-retrieved cloud-top height, 
assuming a cloud thickness of 50 hPa and a positive 
atmospheric lapse rate. Radiosonde data and a forward 
radiative transfer model are used to determine the re- 
lationship between C : and T b : 
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c = — ~ — it . (5) 

‘ iB„ + 5,^ + « 2 W7 - ByW?) 

Tests by Gupta et ah (1992) using products from the 
International Satellite Cloud Climatology Project 
(ISCCP) dataset show that (5) works well except when 
P y - P <h < 200 hPa (/\ is the surface pressure, and P. h 
is the pressure of the cloud base), which is a significant 
problem in the Arctic where low, thin clouds predom- 
inate. 

h. Frouin et at. { 19KH\ 

This algorithm comprises two techniques (hereinafter 
FI and F2) for determining DLF during nighttime. The 
FI method requires temperature profiles, water vapor 
profiles, and cloud properties (cloud-top height and 
cloud fraction), which we obtain from the Path-P dataset 
(Frouin et al. (1988) used operational NOAA TOVS 
retrievals). In this algorithm, clouds are assumed to exist 
in one opaque layer that is 50 hPa thick. The temperature 
and water vapor profiles (from Path-P) are input to 
Streamer to compute the clear-sky DLF To estimate the 
cloudy -sky flux, we use the Path-P cloud effective frac- 
tion and cloud-top height and assign a cloud-base height 
that is 50 hPa lower than the cloud-top height. Using 
the Path-P temperature profile, the cloud-base height is 
matched with the level of the corresponding temperature 
level to obtain the cloud-base temperature. This infor- 
mation is input to Streamer to calculate the cloudy-sky 
DLF. The clear-sky and cloudy-sky fluxes are combined 
according to ( t ). 

In F2, DLF is parameterized as a function of the clear- 
sky flux DLF vlr : 

DLF - DLF lr + o4 , (6) 

which is calculated as for FI . Again, we use A 1>L . in place 
of cloud fraction. The coefficient c depends on latitude, 
season, and cloud type, which is determined by Frouin 
et al. (1988) by simulating DLF in varying cloud con- 
ditions using the Stephens (1978) model. A value of 66 
W m 2 is selected for this study, which corresponds to 
subarctic winter conditions and liquid water clouds. 

Polar clouds rarely exist as a single layer (e.g.. 
Schweiger and Key 1992); thus, we investigate the ef- 
fects of multilayering by applying a cloud-overlap tech- 
nique to the FI method using TOVS Path-P data as 
input. This random-overlap technique is adopted from 
Tian and Curry ( 1989) and is being tested at the NASA 
Langley Research Center for the WCRP-GEWEX SRB 
program. 

The overlap method comprises the following steps: 
) ) obtain cloud-top height and A c , for each satellite re- 
trieval in a 2.5° X 2.5° region; 2) categorize the cloud 
type for each retrieval based on the ISCCP cloud-height 
definitions (high: top pressure below 440 hPa, middle: 
between 440 and 680 hPa. and low: greater than 680 
hPa); 3) determine the fraction of each cloud type in a 


region and F t for high, middle, and low); 4) 

compute the probabilities for each cloud type 

using the following equations: 


C„ = F„n. 

(7) 

C M = |F„/<! - F„)l and 

(8) 

C, = F,H\ - F„ - FJ: 

(9) 


5) calculate clear-sky probabilities (e.g., I - C H for high 
clouds); 6) calculate fractional values for combinations 
of conditions (dear, high alone, high over middle, high 
over low, high over middle over low, etc.); 7) calculate 
fluxes for each case (DLF,) with appropriate cloud-base 
height (50-hPa thickness assumption) and A CL using 
Streamer; and 8) multiply fluxes for each case by their 
corresponding fractional values C,: 

DLF mlp = DLF, C x + DLF 2 C\ + • • • + DLF >V C,. (10) 

where DLF mlp is the new flux value from the cloud- 
overlap method. 

c. Francis (1997) 

The only assumptions in this method are that clouds 
exist in one layer and that cloud fraction is always 100%, 
with all the variability in occurring in the emissivily. 
Clouds may be optically thin and may have varying 
geometric thicknesses. This technique ingests Path-P at- 
mospheric temperature and moisture profiles, effective 
cloud fraction, cloud-top height, and surface tempera- 
ture. Differences between brightness temperatures (TB) 
in several pairs of HIRS channels are used to estimate 
cloud type (positive or negative internal lapse rate), 
phase, thickness, and LWC of Arctic clouds. Cloud 
phase is inferred by subtracting TBs in two pairs of 
channels: HIRS 10 (8.3 gm) and HIRS 8 (11.1 jutm), 
and HIRS 18 (4.0 /tm) and HIRS 19 (3,7 gm). For 
example, in the first pair, the absorption coefficients k uhs 
for water and ice are different. At 8.3 /^m, is similar 
for water and ice, but at 1 1 . 1 gm the difference in 
is large, thereby differentiating ice and water clouds. 

Cloud thickness is estimated using TB differences in 
two pairs of channels: HIRS 6 (13.7 /um) - HIRS 15 
(4.46 /um). and HIRS 14 (4.52 jxm) - HIRS 7 (13.7 
/im). The first pair is for mid- and high clouds (top 
height >750 hPa), and the second pair is used for low 
clouds. Because the weighting function peak of HIRS 
6 is at a lower altitude than that of HIRS 15, its TB is 
warmer in a cloud-free sky. When a thin cloud is present, 
the difference in TB decreases. To determine the cloud 
thickness from the differences in TB, the base fraction, 
a value between 0 and 1, is determined by setting end- 
point thresholds and interpolating linearly between them 
by matching calculated DLFs to observed quantities. 

The liquid or ice water content (IWC) is estimated 
using empirical relationships between mean cloud tem- 
perature and LWC or IWC for water or ice clouds. All 
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this information is input to Streamer to compute DLF. 
See Francis ( 1997) for further details of this algorithm. 

The following three algorithms are simple, empiri- 
cally derived parameteri/ations. For this study the Path- 
P effective cloud fraction is used in place of cloud frac- 
tion in the relationships. 

d. Marshunova ( 1966) 

This method is an empirical y derived parameteri- 
zation to estimate DLF based on surface temperature, 
near-surface vapor pressure, and cloud fraction. A sim- 
ple cloud factor is defined that i icludes the cloud frac- 
tion and a coefficient: 

DLF = DLF 0 ,(I F xA t >, (II) 

where 

DLF clt = oTt(0.67 4 0.0.V”), (12) 

x is a coefficient derived using time-varying surface 
temperatures T A , and e is near^urface vapor pressure. 
In this study x = 0.26 after an analysis by Jacobs ( 1 978). 
Effective cloud fractions and sk n temperatures are ob- 
tained from Path-P data, and v i ; calculated from Path- 
P moisture profiles. 

e. Zi l Itf uui ( 1972) 

This parameterization is a function of both the cloud 
amount and the near-surface air temperature. The Path- 
P surface skin temperatures are used in this study, be- 
cause out analysis of CEAREX measurements reveals 
that the surface skin temperature rarely differs from the 
2-m air temperature by more than 2 K except in pro- 
longed clear winter conditions. 

DLF - DLF tlr + 1^0.96(1 - 9.2 X 10 '•7;)A |, < 13) 

where DLF Ur = oTf (9.2 X 10 *7;). This relationship 
was derived by Zi liman (1972) fiom measurements over 
Antarctic sea ice obtained from Pease (1975). 

f Maykut and Church (1971) 

The relationship was developed with year-round sur- 
face temperature and cloud fra :t ion data collected in 
Barrow. Alaska, over a 5-yr period, during which the 
surface temperature ranged from 244 to 277 K. The DLF 
is parameterized as 

DLF - DLF <)f ( 1 -f 0.22AP), (14) 

where DLF C| , - 0.7855<r7T 

5. Results and discussion 

a. Performance of DLF algor it urns 

Downward longwave fluxes at the surface are com- 
puted using each of the eight rmthods described in sec- 


tion 4 and listed in Table 2. Fluxes calculated from daily- 
average input data are compared to daily-average mea- 
sured fluxes from the CEAREX and LeadEx field pro- 
grams. Scatterplots that illustrate direct comparisons of 
the computed and measured fluxes are shown in Fig. 5, 
and a summary of the comparison statistics is presented 
in Table 3. 

The rms differences and correlation coefficients are 
remarkably similar for all eight methods; hence, they 
are listed in decreasing accuracy according to bias. All 
the methods exhibit negative bias, that is. calculated 
fluxes are too small, although it is negligible for the 
Francis ( 1997) algorithm, which is one of the four spe- 
cifically designed for polar conditions. The cloud-over- 
lap method applied to FI clearly improves the results: 
the bias is reduced from “34 (169F) to -11 (57r) W 
m Results from the other three Arctic-specific meth- 
ods (Marshunova (1966); Maykut and Church (1973); 
Zillman (1972)), although only simple parameteri/a- 
tions. are disappointing. These same three parameteri- 
zations were evaluated by Key et ul. ( 1996) using val- 
idation data from two land stations: Resolute, Northwest 
Territories, Canada, and Barrow. Alaska. When com- 
pared with Key et al. (1996), the Maykut and Church 
(1973) algorithm exhibited a bias one-half as large as 
in our study, probably because it was developed with 
data front Barrow, Alaska. The Zillman (1972) param- 
eterization. developed with data from the Southern 
Ocean was. not surprisingly, the least accurate of the 
three in both evaluations. Our results yielded a larger 
negative bias for the Marshunova (1966) method, again 
probably because the parameterization was developed 
using data from Arctic coastal stations such as Resolute 
and Barrow rather, than observations from within the 
Arctic Ocean. 

h. Sources of error 

1 ) ATMO.SPMF.Rir PROFILES 

The reported bias in Path-P temperature profiles is 
approximately I K (Schweiger et al. 2001 ). which trans- 
lates to an error in DLF of about 3 W m : in summer 
and 1.5 W m : in winter. The bias in water vapor pro- 
files, as compared with radiosonde data from the SHE- 
BA field program, is about 107r. w'hich would produce 
an error in DLF of approximately 3 W in 2 . Wc con- 
sequently conclude that errors in satellite-retrieved tem- 
perature and water vapor profiles do not contribute sig- 
nificantly to the apparent biases in computed DLF. 

2) Cloud fraction 

Because all of these algorithms require information 
about cloud fraction and because DLF in polar regions 
is sensitive to cloud fraction, this variable may account 
for much of the error in computed DLF. As already 
mentioned, the Path-P product used in the analyses is 
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Tabu; 3. Ranked comparison statistic^ of algorithm performance. 


Algorithm reference 

Bias (W m ’) Rinse 
{ A of mean) (W m ’ i 

r* 

Francis (1997) 

-1 (0.1) 

19 

0.77 

Frouin et al. ( 1988). FI with overlap 

II (5) 

19 

0.76 

Frouin et al. (1988). F2 

- 12 (ft) 

19 

0.77 

Marshunova ( 1 966 1 

- 13 (6) 

20 

0.73 

Mavkut and Church (1973) 

- 23 (ID 

20 

0.72 

Gupta (1989) 

- 26 t!2) 

19 

0,75 

Zillman (1972) 

- 30(14) 

22 

0.72 

Frouin et al. (1 988). F! 

34 (I6i 

2(1 

0.73 


* Here r is correlation coefficient. 


the effective cloud traction, which is the product of the 
cloud emissivity and the cloud fraction. Four of the 
algorithms we evaluate assume that clouds are black 
(i.e„ opaque) and 50 hPa thick. These separate as- 
sumptions will contribute to er ors of different signs. 
Clouds with emissivities less tha l unity (nonblack) have 
been observed in the Arctic fCuiry et al. 1996). In these 
cases, which theoretically would yield a retrieved 
that is smaller than the surface-observed cloud fraction, 
the opaque assumption will cause an overestimate in 
DLF because a black cloud will emit more infrared ra- 
diation. all other characteristics being equal. The as- 
sumption of a 50-hPa-thick cloud, on the other hand, 
will cause an underestimation o DLF if the cloud base 
is actually lower. 

To produce a' DLF that is biased by -30 W m \ the 
cloud fraction would have to bt too small by approxi- 
mately 60% for low clouds an J even larger for high 
clouds, according to our sensitivity tests. A comparison 
of cloud climatological descriptions from nine different 
sources ] seven based on human observers and two from 
satellites (Chiacchio 200 1 )| shov's that mean Path-P val- 
ues for April (48%) fall in the middle of the range of 
values [29% (Vowinckel and Or\ ig 1970) to 89% (Barry 
et al. ! 987)1. Excluding (including) the single large out- 
lier, the mean is 46% (51%) anc the standard deviation 
is 8.5% (16.5%). 

In addition, we compare ceilometer-derived cloud 
fractions from the Lead Ex field program to those from 
Path-P (Fig. 6). Although the day-to-day variability is 
lower in Path-P retrievals, the} appear to be slightly 
larger than the ceilometer values. As previously men- 
tioned, however, the ceilometer often did not sense thin 
low clouds, which probably contributes to the large dis- 
crepancies near year days 89, *>9, and 106-109 when 
the ceilometer retrievals were char (Fig. 6c). Data from 
the SHEBA field project analysed by Schweiger et al. 
(2001) did not exhibit this behavior, however; lidar- 
derived cloud fractions were sy Uematicully larger than 
those reported by human observers. Based on these re- 
sults, we conclude that errors n Path-P A ct retrievals 
are not responsible for the large negative biases exhib- 
ited by most of the DLF algorithms. 


3) Cloud-bash hlight 

Cloud-base height in three of the algorithms is de- 
termined by assuming it is 50 hPa lower than the Path- 
P-retrieved cloud-top height. Thus there are two com- 
ponents of this variable to consider: First is the accuracy 
of cloud-top height in the Path-P dataset and how it 
differs from heights estimated using lidar ceilomcters. 
The cloud top can be difficult to define, because cloud 
boundaries are often ephemeral and partially transparent 
to infrared radiation. A comparison of surface-based, 
lidar-radar-retricved cloud-top heights with those from 
Path-P during SHEBA, for example, shows that cloud 
tops are generally higher in lidar-radar retrievals than 
those from Path-P (Schweiger et al. 2001 ). This behav- 
ior is expected because of inherent differences in the 
two observing techniques: the lidar-radar system is sen- 
sitive to the small, sparse ice particles that frequently 
compose high-latitude cloud tops, whereas Path-P re- 
trieves a value corresponding to the effective radiating 
height, that is the height from which the bulk of the 
radiation is emitted from the cloud lop. Although no 
conclusion can be drawn at this time regarding the ve- 
racity of Path-P cloud-top height retrievals, we do know' 
that to contribute to negative biases in computed DLF. 
the retrieved cloud tops would have to be consistently 
too high, which is not what Schw'eiger et al. (2001) 
show. 

The second issue is the assumption of a constant, 50- 
hPa cloud thickness, as in Gupta (1989). FI, FI with 
overlap, and F2. We compare eeilometer-observed 
cloud-base heights from LeadEx to those estimated us- 
ing the Path-P retrieved cloud-top heights assuming a 
50-hPa thickness (Fig. 6a). Cloud-base height observed 
by the ceilometer is markedly lower than that obtained 
assuming a 50-hPa thickness (bias = 1200 m), which 
can be explained either by retrieved cloud tops that are 
too high or thicknesses that are too thin. Whichever the 
cause, this positive bias results in surface fluxes being 
much lower than observed and is the most likely source 
of error in the calculated DLFs. Furthermore, if Path-P 
cloud-top heights arc generally lower than the actual 
values, as suggested by the SHEBA comparison, the 50- 
hPa assumption may be even less realistic than these 
results indicate. Including the cloud-overlap technique 
in FI makes a considerable improvement (Fig. 6b) by 
representing multiple cloud layers and effectively low- 
ering the cloud base. 

Our sensitivity calculations show that if a low-cloud 
base were 50 hPa (about 500 m) too high, the DLF 
would be about 20 W m 2 too small with a typical winter 
Arctic temperature profile. The algorithms that assume 
a 50-hPa-thick cloud exhibit an average bias of ap- 
proximately this amount. We therefore conclude, based 
on these results, that the 50-hPa cloud thickness as- 
sumption is unsuitable for Arctic winter conditions and 
that a more realistic value would be approximately 2 
times as thick. 
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4) Other sources or error 

The poor results exhibited by the Marshunova ( 1 966), 
Zillman (1972), and Maykut and Church ( 1973) param- 
eterization* are somewhat surprising, given that they 
were developed for polar conditions. Because they are 
so simple, it is difficult to ascertain the cause of the 
errors, but it is likely that they arise because the param- 
eterizations were developed using human-observed 
cloud fractions rather than values derived from remote 
sensing instruments. 

Apparent errors in all methods may arise because of 
the inherent differences in comparing daily averaged, 
point-flux measurements with values computed from 
( 100 km) 2 retrieved atmospheric parameters. In addi- 
tion, a negative (clear sky) bias may be introduced be- 
cause TOYS retrieval boxes w ith greater than 90 % ef- 


fective cloud fraction are rejected (Francis 1997). Errors 
may also result from differences in perspectives by sur- 
face observers and satellites for both cloud height and 
cloud fraction. A surface observer has a bottom-up view, 
whereas satellites look down on the cloud top. Surface 
observations of cloud fraction are usually larger than 
those derived from satellites owing to differences in 
view angles and the sky field of view (Schweiger and 
Key 1992; Chiacchio 2001). 

6. Summary and conclusions 

The dominant component of the Arctic surface energy 
budget during almost one-half of the year is longwave 
radiation; however, its spatial and temporal scales of 
variability are not known well, and a reliable, basinwide 
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dataset is not available, A Jong-*erm dataset of surface 
radiation is needed for a variety of applications in cli- 
mate research. A number of methods with varying de- 
grees of complexity exist for estimating DLF from sat- 
ellite and surface measurements and, although surface 
parameterization schemes were intereompared and eval- 
uated thoroughly bv Key et af. (1996), there is little 
information about the performance of satellite 'bused 
techniques in polar nighttime conditions. In this study 
we attempt to evaluate several of these methods and to 
identify reasons for their apparent shortcomings. 

fr> dry polar conditions, cloud properties primarily 
govern DLF. In particular, our sensitivity tests show that 
DLF is most sensitive to uncertainties in cloud fraction. 
LWC and thickness while being, insensitive to droplet 
radius. All of the techniques w.* tested require an es- 
timate of cloud fraction as input, some require temper- 
ature and moisture information, and one ingests satellite 
brightness temperatures directly. Four of the methods 
assume that clouds are a constant 50-hPa thick, three 
contain no thickness information whatsoever, and all 
algorithms but one assume clouds exist in one layer. 
Only one algorithm allows for varying cloud thickness 
and emissivity. Using satellite-derived Path-P products 
as input for the eight DLF methods, we found that all 
techniques except one exhibit t large (>I0 W m ') 
negative bias as compared with measurements from two 
field programs in the Arctic win er and spring. The rms 
errors and correlation eoeffieiei ts do not vary signifi- 
cantly among the methods, however, suggesting that 
cloud fraction, the same data foi which are ingested by 
all methods, probably accounts for much of the vari- 
ability, Our efforts, consequently, focused on identify- 
ing the cause(s) of the consistent negative biases. 

Errors in atmospheric tempeiature and water vapor 
profiles from the Path-P dataset were eliminated as prob- 
able sources of deficient DLLs, Keausc reported biases 
in the Path-P profiles could acccunt for only about 4 W 
m 2 (2% of mean DLF). Path-P t fleet ive cloud fractions 
were also dismissed as a likely source, because com- 
parisons with surface observations (both by humans and 
active remote sensing instruments) reveal a small pos- 
itive bias in cloud fraction, which would result in a 
positive bias in DLF Several of the methods require a 
cloud-top height estimate, for which we used the Path- 
P product. This is a difficult variable to verify owing 
to inherent differences among observing methods, but 
it is likely that Path-P cloud-top heights are biased some- 
what negatively, which would again result in cloud ba- 
ses that were too warm and DLFs that were too large. 
Thus errors in cloud-top height retrievals probably do 
not account for the negatively biased calculated DLF 
values. The most likely source of error, therefore, is the 
assumption by one-half of the algorithms that clouds 
are 50 hPa thick. Model calculat ons show that if a cloud 
layer were 50 hPa too thin, the DLF would be approx- 
imately 20 W m 2 too small. Having eliminated the 
other likely sources of error, w; conclude that the 50- 


hPa-thickness assumption should not he applied during 
polar winter conditions and that, if a constant could 
thickness is required, a value approximately 2 times as 
large should be used. This conclusion is further sup- 
ported by the results of applying a random-eloud-over- 
lap method, which effectively lowers the mean cloud 
base, and by the lack of a DLF bias in the one method 
that ingests the cloud-top height but attempts to estimate 
cloud-base heights from differences in satellite-ob- 
served brightness temperatures. 

Our detailed conclusions are presented below. 

1 ) All eight methods exhibited varying degrees of neg- 
ative bias, ranging from - I (0,5%) to -34 Wm 2 
(16%) as compared with surface measurements of 
DLF; rms errors and correlation coefficients did not 
vary significantly. 

2) The assumption of a constant. 50-hPa cloud thick- 
ness is the most likely cause of negative DLF biases. 
Based on our analysis, we suggest that if a constant 
cloud thickness is used for Arctic w inter conditions, 
it should be doubled to 1 00 hPa, We note, however, 
that the sample size for this analysis is not large (78 
total collocations from CEAREX and LeadEx) and 
validation data cover only one aulumn-w inter season 
in the region northeast of Spitsbergen and one spring 
season in the Beaufort Sea. 

3) Application of a technique to simulate random cloud 
overlapping reduces the bias by the FI method from 
-34 tu “II W m 2 . This improvement is believed 
to be the result of effectively lowering the mean 
cloud base. 

4) The most accurate DLF algorithm may be Francis 
(1997), which exhibits a bias of -1 W m 2 when 
compared with surface measurements. This method 
is specifically developed for Arctic conditions, al- 
lows for clouds within inversion layers, differentiates 
ice and water clouds, and does not assume cloud 
thickness to be constant. We recommend its use for 
Arctic autumn, winter, and spring conditions, its per- 
formance in the summer has not been evaluated yet. 

5) The F2 method outperforms the FI algorithm be- 
cause the coefficient in the formulation depends on 
location, season, and cloud type. 

6) Simple parameteri/.ations designed specifically for 
polar conditions jMaykut and Church (1973); ZilL 
man (1972); Marshunova { 1966)1 did not perform 
well, probably because they were formulated using 
surface-observed cloud fractions from land stations, 
whereas we used satellite-retrieved effective cloud 
fractions, 

7) Errors in satellite-retrieved temperature and water 
vapor profiles do not contribute significantly to the 
negative biases exhibited by calculated DLFs. 

8) Comparisons of Path-P effective cloud fraction with 
ceilometer-retrieved values indicate a slight positive 
bias, which would contribute to a positive bias in 
DLF. We therefore rule out errors in cloud fraction 
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as the source of consistent negative biases in cal- 
culated DLFs. 

9) Attempts to validate Path-P cloud-top heights sug- 
gest they may be too low, which would also con- 
tribute to a positive DLF bias. We therefore dismiss 
this variable as a source of negative DLF biases. 

10) DLF in the winter Arctic is most sensitive to cloud 
fraction, LWC, and cloud thickness; cloud droplet 
size has a negligible effect on DLF 

In summary, we have presented quantitative analyses 
of the sensitivity of downwelling longwave fluxes to 
realistic uncertainties in satellite-derived atmospheric 
parameters in typical Arctic winter conditions. The most 
complex of the algorithms we tested, which includes a 
technique to estimate cloud-base height and emissivity, 
produces DLFs that are closest to measured fluxes in 
these conditions. We found that simpler algorithms that 
assume clouds have a constant thickness and have unit 
emissivity perform poorly in the Arctic, but biases in 
their results are significantly improved by doubling the 
assumed cloud thickness. This analysis should be ex- 
tended to include all seasons, perhaps using measure- 
ments from the SHEBA experiment and more locations 
representative of polar conditions. 
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